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Recent synthesis of fluorinated graphene introduced interesting stable derivatives of graphene. In 
particular, fluorographene (CF), namely fully fluorinated chair conformation, is found to display 
crucial features, such as high mechanical strength, charged surfaces, local magnetic moments due 
to vacancy defects and a wide band gap rapidly reducing with uniform strain. These properties, 
as well as structural parameters and electronic densities of states are found to scale with fluorine 
coverage. However, most of the experimental data reported to date neither for CF, nor for other 
CnF structures complies with the results obtained from flrst-principles calculations. In this study, 
we attempt to clarify the sources of disagreements. 

PACS numbers: 73.22.Pr, 63.22.Rc, 61.48.Gh, 71.20.-b 



I. INTRODUCTION 

Active research on graphen^ revealed not only nu- 
merous exceptional propertie^^HD b^t also have prepared 
a ground for the discovery of several graphene based 
materials. Preparation of freestanding graphene sheets 
with nonuniform oxygen coverage have been achieved.'^ 
More recently the synthesis of two-dimensional hydro- 
carbon in honeycomb structure, so called graphan^ 
(CH). showing diverse electronic, magnetic and mechan- 
ical properties^ is reported. 

According to Pauling scale F has electronegativity of 
3.98, which is higher than those of C(2.55), H(2.20) and 
0(3.44), and hence fluorination of graphene is expected 
to result in a material, which may be even more inter- 
esting than both graphene oxide and CH. Before the 
first synthesis of gra phene , fluorinated graphite has been 
treated theoretically.l^'^ Owing to promising properties 
revealed for CH, fluorinated graphene structures are now 
attracting considerable interesi!i3l22l despite uncertainties 
in their chemical compositions and atomic structures. In 
an effort to identify the structures of fluorinated sam- 
ples, previous theoretical models attempted to deduce 
the lowest energy structures.E^E^ In addition, band gaps 
of different structures calculated within Density Func- 
tional Theory (DFT) are compared with the values re- 
vealed through specific measurements.'^^'iS However, nei- 
ther the stability of proposed structures are questioned, 
nor underestimation of band gaps within DFT has been 
a subject matter. Raman spectrum by itself, has been 
limited in specifying C„F structures. -"^^ 

In this work, we first determined stable C„F struc- 
tures for n < 4. Then we revealed specific proper- 
ties (such as internal structural parameters, elastic con- 
stants, formation and binding energies, energy band gap 
and photoelectric threshold) for those stable structures 
as signatures to identify the derivatives probed experi- 
mentally. We placed an emphasis on fully fiuorinated 
graphene or fluorographene (CF), in which D and G Ra- 
man peaks of bare graphene disappear after long fluori- 
nation period. li^IIl Present study reveals that the proper- 



ties, such as structural parameters, binding energy, band 
gap and phonon modes of various fluorinated structures 
are strongly dependent on the binding structure of F 
atoms and their composition. Some of these properties 
are found to roughly scale with F coverage. While the 
stable C2F chair structure is metallic, CF is a nonmag- 
netic insulator with a band gap. Eg, being much larger 
than 3 eV, i.e. a value attributed experimentally to fully 
fluorinated graphene. In view of the calculated diffu- 
sion constant, Raman active modes and other proper- 
ties, available experimental data suggest that domains 
(or grains) of various C„F structures with extended and 
imperfect grain boundaries can coexist after fluorination 
process. Hence the measured properties are averaged 
from diverse perfect and imperfect regions. 

II. COMPUTATIONAL METHODOLOGY 

Our predictions are obtained from flrst-principles plane 
wave calculations^^ within DFT, which is demonstrated 
to yield rather accurate results for carbon based mate- 
rials. Calculations are performed using spin-polarized 
local density approximation (LDA)'^ and projector aug- 
mented wave (PAW) potentials. Kinetic energy cutoff 
h?\\i + Gp/2TO for plane-wave basis set is taken as 500 
eV. In the self-consistent potential and total energy cal- 
culations of fluorographene a set of (25x25x1) k-point 
sampling is used for Brillouin zone (BZ) integration. 
The convergence criterion of self consistent calculations 
for ionic relaxations is 10~^ eV between two consecu- 
tive steps. By using the conjugate gradient method, 
all atomic positions and unit cells are optimized until 
the atomic forces are less than 0.03 eV/A. Pressures on 
the lattice unit cell are decreased to values less than 0.5 
kBar. The energy band gap, which is usually underesti- 
mated in DFT, is corrected by frequency-dependent GWq 
calculations.^^ In GWq corrections screened Coulomb po- 
tential, W, is kept fixed to initial DFT value Wq and 
Green's function, G, is iterated four times. Various tests 
are performed regarding vacuum spacing, kinetic energy 
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cut-oflF energy, number of bands, k-points and grid points. 
Finally, the band gap of CF is found 7.49 eV after GWq 
correction, which is carried out by using (12x12x1) k- 
points in BZ, 15 A vacuum spacing, default cut-off po- 
tential for GWo, 192 bands and 64 grid points. Phonon 
frequencies and phonon eigenvectors are calculated using 
the Density Functional Perturbation Theory (DFPT).'23 



III. STRUCTURES OF FLUORINATED 
GRAPHENE 

Each carbon atom of graphene can bind only one F 
atom and through coverage (or decoration) of one or two 
sides of graphene, one can achieve diverse C„F structures. 
Uniform F coverage is specified by = 1/n (namely one 
F adatom per n C atoms), whereby Q = 0.5 corresponds 
to half fluorination and 9 = 1 is fluorographene CF. The 
adsorption of a single F atom to graphene is precursor for 
fluorination. When placed at diverse sites of a (4 x 4) su- 
percell of graphene, simple F atom moves to the top site 
of a carbon atom and remains adsorbed there. For the 
resulting structure spin-polarized state (with 0.4 Bohr 
magneton) is only 2 meV favorable than the nonmag- 
netic state. The binding energy of F is Eb=2.1\ eV in 
equilibrium, that is a rather strong binding unlike many 
other adatoms adsorbed to graphene. An energy bar- 
rier, Qb =^0.45 eV, occurs along its minimum energy 
migration path. Our calculations related with the mini- 
mum energy path of a single F atom follow hexagons of 
underlying graphene. Namely, F atom migrates from the 
highest binding energy site, i.e. top site (on top of carbon 
atom) to the next top site through bridge site (bridge po- 
sition between two adjacent carbon atoms of graphene). 
The corresponding diffusion constant for a single F atom, 
D = j/ae~'5^/'^^"^ is calculated in terms of the lattice 
constant, a = 2.55 A and characteristic jump frequency 
V w39 THz. Experiments present evidences that energy 
barrier s on t he order of 0.5 eV would make the adatoms 
mobile .'i^'^ Moreover, this energy barrier is further low- 
ered even it is collapsed in the presence of a second F 
atom at the close proximity. Consequently, this situation 
together with the tendency towards clustering favors that 
C„F grains (or domains) of different n on graphene can 
form in the course of fluorination. We note that the en- 
ergy barrier for the diffusion of a single carbon adatom 
adsorbed on the bridge sites of graphene was calculated 
to be in the similar energy range. Carbon adatoms on 
graphene were found to be rather mobile. That energy 
barrier for single C adatom was found to decrease, even 
to collapse at the close proximity of a second adatom. ■^^ 

In earlier theoretical studies ,'^^'^^^3 the total energies 
and/or binding energies were taken as criteria for whether 
a given C„F structure exists. Even if a C„F structure 
seems to be in a minimum on the Born- Oppenheimer 
surface, its stability is meticulously examined by calcu- 
lating frequencies of all phonon modes in BZ. Here we 
calculated phonon dispersions of most of optimized C„F 
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FIG. 1: (Color online) Atomic structure and calculated 
phonon bands (i.e. phonon frequencies versus wave vector, 
k) of various optimized C^F structures calculated along the 
symmetry directions of BZ. Carbon and fluorine atoms are 
indicated by black (dark) and blue (light) balls, respectively, 
(a) C2F Chair structure, (b) C2F Boat structure, (c) C4F 
structure. Units are A for structural parameters and cm"'^ 
for frequencies. 



structures. We found C4F, C2F boat, C2F chair (See 
Fig. [ij and CF chair (See Fig. [2]) structures have posi- 
tive frequencies throughout the Brillouin zone indicating 
their stability. 

Some of phonon branches of C„F structures (for ex- 
ample, CF boat) have imaginary frequencies and hence 
are unstable, in spite of the fact that their structures can 
be optimized. The possibility that these unstable struc- 
tures can occur at flnite and small sizes is, however, not 
excluded. For stable structures, the gap between opti- 
cal and acoustical branches is collapsed, since the optical 
branches associated with the modes of C-F bonds oc- 
cur at lower frequencies. This situation is in contrast 
with the phonon spectrum of graphane,!^ where optical 
modes related with C-H bonds appear above the acous- 
tical branches at ~ 2900 cm~^. 

The formation energy of fluorination is defined as Ef ^ 
{np^ET,F2 + Et^gt — Et^^f) /np in terms of the total 
ground state energies of optimized structures of graphene 
and fluorinated graphenes at different composition, re- 
spectively, Ex^Gri Et,c„Fi and the total ground state 
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TABLE I: Comparison of the calculated properties of four stable, fluorinated graphene structures (namely CF, C2F chair, C2F 
boat and C4F) with those of graphene and CH. Lattice constant, a — b {a ^ b for rectangular lattice); C-C bond distance, 
dec (second entries with slash differ from the previous one); C-X bond distance (X indicating H (F) atom for CH (CF)), 
dcx; the buckling, 5; angle between adjacent C-C bonds, ac; angle between adjacent C-X and C-C bonds, ax', total energy 
per cell comprising 8 carbon atoms Et; formation energy per X atom relative to graphene, Ef, binding energy per X atom 
relative to graphene, Eb (the value in parenthesis E^^' excludes the X-X coupling); desorption energy, Ed (see the text for formal 
definitions); energy band gap calculated by LDA, Eg^^; energy band gap corrected by GWo, E^^"; photoelectric threshold, 
$; in-plane stiffness, C; Poison ratio, i/. All materials are treated in hexagonal lattice, except C2F boat having rectangular 
lattice. 
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5.99 
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FIG. 2: (Color online) (a) Atomic structure of fluorographene 
CF. a and b are the lattice vectors (|a| = |6|) of hexago- 
nal structure; dec {dcp) is the C-C (C-F) bond distance; 5 
is the buckling, (b) Phonon frequencies versus wave vector 
k of optimized CF calculated along symmetry directions in 
BZ. (c) Symmetries, frequencies and descriptions of Raman 
active modes of CF. (d) Calculated Raman active modes of 
graphene, CH, CF and C4F are indicated on the frequency 
axis. Those modes indicated by " + " are observed experimen- 
tally. There is no experimental Raman data in the shaded 
regions. Units are A for structural parameters and cm~^ for 
frequencies. 



energy of single carbon atom, Et.c, of F2 molecule and 
F atom, Et.F2 Et,f- Similarly, the binding energy 
of F atom relative to graphene including F-F coupling is 
Eh — {Et^gt + npET^F — ETfinP) /np and without F-F 



coupling E^,' = {Et^gt + Et^ufF - ET,G^F)/nF- Here 
Ex^npF is the total energy of suspended single or double 
layers of F occupying the same positions as in C„F. The 
desorption energy, Ed is the energy required to remove 
one single F atom from the surface of C„F. nF2 and uf 
are numbers of F2 molecules and F atoms, respectively. 
The total energies are calculated in periodically repeat- 
ing supercells comprising 8 carbon atoms and keeping 
all the parameters of calculations described above using 
spin-polarized as well as spin unpolarized LDA. Lowest 
(magnetic or nonmagnetic) total energy is used as ground 
state total energy. 

Fluorographene (CF), where F atoms are bound to 
each C atom of graphene alternatingly from top and bot- 
tom sides is energetically most favorable structure. Upon 
full fluorination, the planar honeycomb structure of C 
atoms becomes buckled (puckered) and C-C bond length 
increases by ~10%. At the end, while planar sp^-bonding 
of graphene is dehybridized, the buckled configuration is 
maintained by sp'^-like rehybridization. In Table |l] the 
calculated lattice constants, internal structural parame- 
ters, relevant binding energies and energy band gaps of 
stable C„F structures are compared with those of bare 
graphene and CH.^'' Notably, internal parameters (such 
as (S, C-C bond length) as well as lattice constants of 
various C„F structures vary with F coverage, Q. CF, 
has highest values for Ej, Ei,, E^^i and Ed given in Ta- 
ble |l] those of C4F are second highest among stable C„F 
structures. 

Since Raman spectrum can convey information for a 
particular structure and hence can set its signature, cal- 
culated Raman active modes of stable C4F and CF struc- 
tures together with those of graphene and CH are also 
indicated in Fig. |2] (c) and (d). It is known that the 
only characteristic Raman active mode of graphene at 
1594 cm~^ is observed so far.l^ Similarly, for CH the 
mode at ^1342 cm~^ is observed.'^ One of two Raman ac- 
tive modes of C4F at 1645 cm~^ seems to be observed.!^ 
In compliance with the theory, '^^ phonon branches of all 
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these observed modes exhibit a kink structure. However, 
none of the Raman active modes of CF revealed in Fig. [2] 
has been observed yet. Raman spectroscopy in the low 
frequency range may be useful in identifying experimen- 
tal structures. 



IV. ELECTRONIC STRUCTURES 

Energy bands, which are calculated for the optimized 
C4F, C2F boat, C2F chair and CF chair structures are 
presented in Fig. [3] and Fig. |4] structures, respectively. 
The orbital projected densities of states (PDOS) together 
with the total densities of states of these optimized struc- 
tures are also presented. Analysis of the electronic struc- 
ture can also provide data to reveal the observed struc- 
ture of fluorinated graphene. As seen in Table |T] stable 
C„F structures have LDA band gaps ranging from eV 
to 2.96 eV. Surprisingly, C2F chair structure is found to 
be a metal owing to the odd number of valance electrons 
in the primitive unit cell. Even if various measurements 
on the band gap of fluorinated graphene lie in the en- 
ergy range from 68 meVG^l to 3 eV^ these calculated 
band gaps are underestimated by LDA. Incidentally, the 
band gaps change significantly after they are corrected 
by various self-energy methods. In fact, the correction 
using GWo self-energy method predicts a rather wide 
band gap of 7.49 eV for CF. The corrected band gaps 
for C2F boat structure and C4F are 5.68 eV and 5.99 eV, 
respectively. It should be noted that the GWq self-energy 
method has been successful in predicting the band gaps 
of 3D semiconductors.132 

While predicting much larger band gap for CF, the 
measured band gap of ~3 eV reported by Nair et al.l^ 
marks the serious discrepancy between theory and ex- 
periment. The character of the band structure of CF are 
revealed from the analysis of projected density of states 
as well as charge densities of specific bands in Fig.|4](b). 
Conduction band edge consists of the antibonding com- 
bination of p^-orbitals of F and C atoms. The p^-orbitals 
of C atoms by themselves, are combined to form 7r-bands. 
The bands at the edge of the valence band are derived 
from the combination of C-px +Py and F-p^ ~^Py orbitals. 
The total contribution of C orbitals to the valence band 
can be viewed as the contribution of four tetrahedrally 
coordinated sp'^-like hybrid orbitals of s- and p-orbitals 
of C atoms. However, the deviation from tetrahedral co- 
ordination increases when n increases or single side is flu- 
orinated. As a matter of fact, the total density of states 
presented in Fig. |3] and Fig. |4] mark crucial differences. 
In this respect, spectroscopy data is expected to yield 
significant information regarding the observed structures 
of fluorinated graphenes. 

The contour plots of the total charge density, px, in the 
F-C-C-F plane suggests the formation of strong covalent 
C-C bonds from the bonding combination of two C-sp'^ 
hybrid orbitals. The difference charge density, Ap (which 
is obtained by subtracting the charges of free C and free 
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FIG. 3: (Color online) Energy band structures of various sta- 
ble CnF structures together with the orbital projected den- 
sities of states and the total densities of states (DOS). The 
LDA band gaps are shaded and the zero of energy is set to 
the Fermi level Ep. Total DOS is scaled to 45%. Valence and 
conduction band edges after GWo correction are indicated 
by filled/red circles, (a) C2F chair structure, (b) C2F Boat 
structure, (c) C4F structure. 



F atoms situated at their respective positions in CF), in- 
dicates charge transfer to the middle of C-C bond and to 
F atom, revealing the bond charge between C atoms and 
ionic character of C-F bond. However, the value of charge 
transfer is not unique, but diversifies among different 
methods of analysis.^^ Nevertheless, the direction of cal- 
culated charge transfer is in compliance with the Pauling 
ionicity scale and is corroborated by calculated Born ef- 
fective charges, which have in-plane (||) and out-of-plane 
(_L) components on C atoms, y = 0.30, Z"^ j_ = 0.35 
and on F atoms p = —0.30, Zpj_ = —0.35. 

Finally, we note that perfect CF is a nonmagnetic insu- 
lator. However, a single isolated F vacancy attains a net 
magnetic moment of 1 Bohr magneton (hb) and localized 
defect states in the band gap. Creation of an unpaired tt- 
electron upon F vacancy is the source of magnetic mo- 
ment. However, the exchange interaction between two 
F-vacancies calculated in a (7x7x1) supercell is found to 
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FIG. 4: (Color online) (a) Energy band structure of CF to- 
gether with the orbital projected DOS and total densities of 
states. The LDA band gap is shaded and the zero of energy 
is set to the Fermi level, Ep. Valence and conduction band 
edges after GWo correction are indicated by filled/red cir- 
cles, (b) Isosurfaces of charge densities of states corresponds 
to first (VI), second (V2) valance and first (CI) and second 
(C2) conduction bands at the F- and AT-points. (c) Contour 
plots of the total charge density pr and difference charge den- 
sity Ap in the plane passing through F-C-C-F atoms. Contour 
spacings are 0.03 e/A^. 

be nonmagnetic for the first nearest neighbor distances 
due to spin pairings. Similar to graphane,^*' it is also 
possible to attain large magnetic moments on F-vacant 
domains in CF structures. 



V. ELASTIC PROPERTIES OF CF 

Having analyzed the stability of various C„F structures 
with n =1,2 and 4, we next investigate their mechani- 
cal properties. The elastic properties of this structure 
can conveniently be characterized by its Young's modu- 
lus and Poison's ratio. However, the in-plane stiffness C 



is known to be a better measure of the strength of single 
layer honeycomb structures, since the thickness of the 
layer h cannot be defined unambiguously. Defining Aq 
as the equilibrium area of a C„F structures, the in-plane 
stiffness is obtained as C = {6E'^/Se^)/Ao, in terms of 
strain energy Eg and uniaxial strain e.^^ The values of 
in-plane stiffness C, and Poisson's ratio v, calculated for 
stable C„F structures are given in Table |T] together with 
the values calculated for graphene and graphane. For 
example, the calculated values of CF are C=250 J/m^ 
and i/=0.14:. It is noted that C increases with n. For 
CF (i.e. n = 1), the in-plane stiffness is close to that 
calculated for CH. It appears that either the interaction 
between C-F bonds in CF (or the interaction between 
C-H bonds in CH) does not have significant contribution 
to the in-plane stiffness. The main effect occurs through 
dehybridization of sp^ bonds of graphene through the 
formation C-F bonds (or C-H bonds). 

A value of Young's modulus around 0.77 TPa can be 
calculated by estimating the thickness of CF as /i=3.84 
A, namely the sum of the thickness of graphene (3.35 A) 
and buckling, S (0.49 A). This value is smaller, but com- 
parable with the value proposed for graphene, i.e. ^ 1 
TPa. Here the contribution of C-F bonds to the thick- 
ness of CF is neglected, since the interaction between C-F 
bonds has only negligible effects on the strength of CF. 

In Fig. |5] the variation of strain energy Eg and its 
derivative, SE^/Se, with strain, e are presented in both 
elastic and plastic regions. Two critical strain values, 
and ec2, are deduced. The first one, e^, is the point 
where the derivative curve attains its maximum value. 
This means that the structure can be expanded under 
smaller tension for higher values of strain. This point 
also corresponds to phonon instabilitj^^l where the lon- 
gitudinal acoustic modes start to become imaginary for 
e > Eci- The second critical point, (~ 0.29), corre- 
sponds to the yielding point. Until this point the honey- 
comb like structure is preserved, but beyond it the plastic 
deformation sets in. We note that for < e < the 
system is actually in a metastable state, where the plastic 
deformation is delayed. Under the long wavelength per- 
turbations, vacancy defects and high ambient tempera- 
ture approaches to Cc^. In fact, our further molecular 
dynamics simulations show that 0.17 at 300K and 

to 0.16 at 600K. In the presence of periodically repeat- 
ing F-vacancy and C+F-divacancy, the value of is 
also lowered to 0.21 and 0.14, respectively. Apart from 
phonon instability occurring at high e, the band gap is 
strongly affected under uniform expansion. In Fig. 
we show the variation of LDA and GWq- corrected band 
gaps under uniform expansion. The LDA gap slightly 
increases until e = 0.05 and then decreases steadily with 
increasing e. The GW^o-corrected band gap essentially 
decreases with increasing strain. For example, its value 
decreases by 38% for e = 0.20. 
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Strain e [Ac/c^j] Strain e [Ac/c^] 

FIG. 5: (Color online) (a) Variation of strain energy and its first derivative with respect to the uniform strain e. Orange/gray 
shaded region indicates the plastic range. Two critical strains in the elastic range are labeled as tci and ec2. (b) Variation of 
the band gaps with e. LDA and GWo calculations are carried out using 5x5 supercell having the lattice parameter of co=5a, 
and Ac is its stretching 



VI. CONCLUSIONS 

Present analysis of fluorinated graphenes shows that 
different C„F structures can form at different level of F 
coverage. Calculated properties of these structures, such 
as lattice parameter, dec distance, band gap, density of 
states, work function, in plane stiffness C, Poisson's ratio 
and surface charge, are shown to depend on n or coverage 
Q. Relevant data reported in various experiments do not 
appear to agree with the properties calculated any one 
of the stable C„F structures. This finding lets us to con- 
clude that domains of various C„F structures can form 
in the course of the fluorination of graphene. Therefore, 
the experimental data may reflect a weighted average of 
diverse C„F structures together with extended defects in 
grain boundaries. In this respect, imaging of fluorinated 
graphene surfaces by scanning tunneling and atomic force 
microscopy, as well as x-ray photoemission spectroscopy 
is expected to shed light on the puzzling inconsistency 
between theory and experiment. 

Finally, our results show a wide range of interesting 
features of C„F structures. For example, a perfect CF 
structure as described in Fig. [2] is a stiff, nonmagnetic 



wide band gap nanomaterial having substantial surface 
charge, but attains signiflcant local magnetic moment 
through F-vacancy defects. Moreover, unlike graphane, 
half fluorinated graphene with only one side fluorinated 
is found to be stable, which can be further functional- 
ized by the adsorption of adatoms to other side. For 
example, hydrogen atoms adsorbed to other side attain 
positive charge and hence permanent transversal electric 
fleld, which can be utilized to engineer electronic proper- 
ties. 

VII. ACKNOWLEDGMENTS 



This work is supported by TUBITAK through Grant 
No:108T234. Part of the computational resources has 
been provided by UYBHM at ITU through grant Grant 
No. 2-024-2007. We thank the DEISA Consortium 
(www.deisa.eu), funded through the EU FP7 project RI- 
222919, for support within the DEISA Extreme Com- 
puting Initiative. S. C. acknowledges the partial support 
of TUBA, Academy of Science of Turkey. The authors 
would also hke to acknowledge the valuable suggestions 
made by D. Alfe. 



* Electronic address: ciraci@fen.bilkent.edu.tr 

^ K. S. Novoselov, A. K. Geim, S. V. Morozov,D. Jiang, 

Y. Zhang, S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, 

Science 306, 666 (2004). 
^ A. K. Geim, K. S. Novoselov, Nat. Mater. 6, 183 (2007). 
^ C. Berger, Z. Song, T. Li, X. Li, A. Y.Ogbazghi, R. Feng, 

Z. Dai, A. N. Marchenkov, E. H. Conrad, P. N. First, and 

W. A. de Heer, Science 312, 1191 (2006). 

* M. L Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. 
Phys. 2, 620 (2006). 

^ H. §ahin, R. T. Senger and S. Ciraci, J. Appl. Phys. 108, 
074301 (2010). 



" D. A. Dikin, S. Stankovich, E. J. Zimney, R. D. Finer, G. 
H. B. Dommett, G. Evmenenko, S. T. Nguyen and R. S. 
Ruoff, Nature 448, 457 (2007). 

D. C. Ellas, R. R. Nair, T. M. G. Mohiuddin, S. V. 

Morozov, P. Blake, M. P. Halsall, A. C. Ferrari, D. W. 

Boukhvalov, M. L Katsnelson, A. K. Geim and K. S. 

Novoselov, Science 323, 610 (2009). 
^ J. O. Sofo, A. S. Chaudhari and G. D. Barber, Phys. Rev. 

B 75, 153401 (2007). 
® D. W. Boukhvalov, M. L Katsnelson, A. L Lichtenstein, 

Phys. Rev. B 77, 035427 (2008). 
" H. §ahin, C. Ataca and S. Ciraci, Appl. Phys. Lett. 95, 



7 



222510 (2009). 

" H. §ahin, C. Ataca and S. Ciraci, Phys. Rev. B 81, 205417 
(2010). 

M. Topsakal, S. Cahangirov and S. Ciraci, Appl. Phys. 
Lett. 96, 091912 (2010). 

J.-C. Charlier, X. Gonze and J. -P. Michenaud, Phys. Rev. 
B 47, 16152 (1993). 

Y. Takagi and K. Kusakabe, Phys. Rev. B, 65, 121103 

(2002). 

1^ D. W. Boukhvalov, Physica E, 43, 199 (2010). 

1® S. -H. Cheng, K. Zou, F. Okino, H. R. Gutierrez, A. Gupta, 

N. Shen, P. C. Eklund, J. O. Sofo, and J. Zhu, Phys. Rev. 

B 81, 205435 (2010). 
^"^ J. T. Robinson, J. S. Burgess, C. E. Junkermcier, S. C. 

Badescu, T. L. Reinccko, F. K. Perkins, M. K. Zalalutd- 

niov, ,J. W. Baldwin, J. C. Culbertson, P. E. Sheehan, and 

E. S. Snow, Nano Lett. 10, 3001 (2010). 

R. R. Nair, W. Ren, R. Jahl, L Riaz, V. G. Kravets, L. 

Britnell, P. Blake, F. Schedin, A. S. Mayorov, S. Yuan, M. 

I. Katsnelson, H.-M. Cheng, W. Strupinski, L. G. Bulu- 

sheva, A. V. Okotrub, I. V. Grigorieva, A. N. Grigorenko, 

K. S. Novoselov, A. K. Geim, Small 6, 2877 (2010). 

0. Leenaerts, H. Peelaers, A. D. Hernandez-Nieves, B. Par- 
toens, and F. M. Peeters, Phys. Rev. B 82, 195436 (2010). 

^° F. Withers, M. Dubois, and A. K. Savchenko, Phys. Rev. 
B 82, 073403 (2010). 

M. Klintenberg, S. Lebegue, M. I. Katsnelson, and O. 
Eriksson, Phys. Rev. B 81, 085433 (2010). 

22 E. Munoz, A. K. Singh, M. A. Ribas, E. S. Pencv, and B. 

1. Yakobson, Diamond Relat. Mater. 19, 368 (2010). 

23 G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993); G. 
Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 (1996). 

2* D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566, 
(1980). 

25 P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
2® M. Shishkin and G. Kresse, Phys. Rev. B 74, 035101 
(2006). 



P. Giannozzi et al. J. Phys. Condens. Matter 21, 395502 
(2009). 

2« Y. Can, L. Sun, and F. Banhart, Small 4, 587 (2008). 
2^^ C. Ataca, E. Akturk, H. §ahin and S. Ciraci, J. Appl. Phys. 

109, 013704 (2011). 
3° A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, M. 

Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. Novoselov, 

S. Roth, and A. K. Geim, Phys. Rev. Lett. 97, 187401 

(2006). 

3^ S. Piscanec, M. Lazzeri, F. Mauri, A. C. Ferrari, and J. 
Robertson, Phys. Rev. Lett. 93, 185503 (2004). 

32 HSE hybrid functional implemented in VASP [K. Hum- 
mer, J. Harl, G. Kresse, Phys. Rev. B 80, 115205 (2009)] 
is demonstrated to be as succesful as GWo self-energy cor- 
rection method in predicting the band gaps of bulk 3D 
crystals. Consistently, for 2D honeycomb structures HSE 
is found to yield smaller values than those of GWq. For 
example, while HSE predicts the band gap of CF as 4.86 
eV, GWo gives 7.49 eV. Similar trends are also found for 
CH and 2D Boron Nitride (BN). HSE and GWo corrected 
band gaps of CH (BN) band gaps are 4.51 (5.74) and 5.97 
(6.86) eV, respectively. 

33 H. §ahin, S. Cahangirov, M. Topsakal, E. Bekaroglu, E. 
Akturk, R. T. Senger, and S. Ciraci, Phys. Rev. B 80, 
155453 (2009). 

3^ For example, Bader [G. Henkelman, A. Arnaldsson, H. 
Jonsson, Comput. Mater. Sci.36, 254-360 (2006)], Lowdin 
JR-O. Lowdin, J. Chem. Phys. 18, 365 (1950)] and Mul- 
liken [R. S. MuUiken, J. Chem. Phys. 23, 1841 (1955)] 
methods predict charge transfer from C to F atoms, re- 
spectively, 0.59, 0.19 and 0.06 electrons. Further analysis 
by calculating planarly averaged charge density of CF and 
of free C and F atoms, where the charge due to the tails 
of orbitals are carefully accounted deduces that the total 
of 0.11 electrons are transferred from buckled graphene to 
each F atom. 



